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We study the time-dependent transport of charge and spin through a ring-shaped region sequen- 
tially coupled to a weakly interacting quantum dot in the presence of an Aharonov-Bohm flux 
and spin-orbit interaction. The time-dependent modulation of the spin-orbit interaction, or of the 
corresponding Aharonov-Casher flux, together with the modulation of the dot-level induces an elec- 
trically pumped spin current even in absence of a charge current. The results beyond the adiabatic 
regime show that an additional rectification current proportional to cos(y>), being ip the relative 
phase between the time varying parameters, is generated. We discuss the relevance of such term in 
connection with recent experiments on out-of-equilibrium quantum dots. 



I. INTRODUCTION 

In recent years the process of miniaturization of man- 
made electronic circuits has permitted to reach the 
molecular scale providing the opportunity of testing 
quantum mechanics in nano-electronics measurements^. 
In this framework, one of the most exciting challenges 
is to encode information by means of the electron spin 
instead of the charge, giving rise to the so-called spin- 
based electronics or spintronics^,. To face the need of 
spin-polarizing systems acting as a source in spintron- 
ics devices, one promising possibility is to exploit the 
quantum interference effects by external electric or mag- 
netic fields. In ring-shaped structures made of semicon- 
ducting materials a spin sensitive phase, the Aharonov- 
Casher phased, is originated by the Rashba spin-orbit 
interaction^. Such phase combined with the Aharonov- 
Bohm phase 7 induced by a magnetic field is an useful 
tool to achieve spin-polarizing devices^. Another in- 
teresting phase interference effect is originated by the 
periodic modulation of two out-of-phase parameters af- 
fecting the scattering properties of a nanostructure. Such 
phase effect, known in the literature as quantum pump- 
ing, was first introduced by Thoulessi£. After the Thou- 
less theory, a scattering approach to the adiabatic quan- 
tum pumping was formulated by P. W. Brouwer— who 
showed that the d.c. current pumped by means of an 
adiabatic modulation of two out-of-phase independent 
parameters can be expressed in terms of the paramet- 
ric derivatives of the scattering matrix. In the adiabatic 
regime described by the Brouwer formula, i.e. when the 
pumping frequency is much slower than the tunneling 
rates, a d.c. current proportional to u> sin(</?) is orig- 
inated, being ip the phase difference between the two 
parameters. Such theoretical prediction has been veri- 
fied experimentally by M. Switkes et alJ^ even though 
some anomalies in the current-phase relation have been 
reported. In particular, it has been observed a non- 
vanishing current at ip = 0. Several anomalies observed 
in the experiment can be explained by rectification of 
a.c. displacement currents as proposed in Ref.[l3j]. Ac- 



cording to this work, the rectified currents are responsi- 
ble for measurable effects which may be dominant over 
the pumping currents. In order to discriminate between 
rectified currents and pumping effects symmetry argu- 
ments can be exploited. For instance, the d.c. rectifi- 
cation voltage V rec t is symmetric under reversal of the 
magnetic field V rec t{B) = V rec t(—B), while the voltage 
generated by a quantum pum p is not. On the other hand, 
it has been shown in Ref.[14[ that finite frequency effects, 
considered within a non-equilibrium Green's function ap- 
proach, can lead to current-phase relations of the form 
I c ~ a shi(ip)-\-bcos((p), where the coefficients a and b are 
function of the pumping frequency to. Differently from 
the adiabatic regime where the pumped currents are odd 
function of the relative phase <p> between the modulated 
parameters, any symmetry can be realized in the general 
non-equilibrium case. 

In the following we analyze the charge and spin pump- 
ing in a ring-shaped conductor sequentially coupled to a 
quantum dot (see Fig[l} and apply the non-equilibrium 
Green's function approach to analyze the dc current from 
the adiabatic to non-adiabatic regime addressing the 
question about the existence of rectification terms. In the 
ring region shown in Fig[T]the electrons feel an Aharonov- 
Bohm phase associated to a time varying Aharonov- 
Casher phase. The last is related to the Rashba spin- 
orbit interaction which is tunable by means of a gate 
voltage 15 . An additional time-dependent modulation of 
the dot energy level is also considered . If no voltage bias 
is present between the two external leads, the electron 
current is activated by absorption and emission of quan- 
tized photon energy. Thus, in the following the charge 
and spin pumped current are studied as a side-effect of 
boson-assisted tunneling. 

The paper is organized as follows. In section [TT] we intro- 
duce the model Hamiltonian and derive the general ex- 
pression for the non-equilibrium Green's function and re- 
spective self-energies for the non-interacting and weakly 
interacting case. In section Hill we employ a one-photon- 
approximation and obtain a compact expression for the 
d.c. current pumped in the left lead. In section IIV1 we 



2 



□ 




Leads 

Interferometer 
|^ Quantum Dot 

FIG. 1: An Aharonov-Bohm-Casher quantum ring sequen- 
tially coupled to a quantum dot. The energy on the dot and 
the Aharonov-Casher flux are modulated in time with fre- 
quency UJ. 



II. THE MODEL AND NON-EQUILIBRIUM 
CURRENT 



The Hamiltonian of an Aharonov-Bohm-Casher ring 
sequentially coupled to an interacting quantum dot in 
the presence of time varying parameters can be written, 
in the local spin frame, as follows^: 

H(t) = H c + Y^ e{t)dld a + Un^ [ n i + 

+ ^[2ucos($ CT (t))4c*; CT i +wcl <Tr d (T ] +H.c,(l) 



where H c = J2k<j, a ={i,r} £ k c kaa c k<ya is the free electrons 
Hamiltonian describing the left/right (l/r) leads kept at 
the same chemical potential /i. The second and third 
term represent the dot Hamiltonian consisting of the 
electron-electron interaction term l/n^ni (n a = d^da) 
and of the time dependent dot energy level e(t) — e + 
sm(u)t + </?), being ui the frequency of the modulation. 
The last term in Hamiltonian describes the tunneling be- 
tween the left lead and the dot, ucos($ cr ), and the right 
lead and the dot through a tunnel barrier, w. The trans- 
mission coefficients u and w, which in general may be spin 
and momentum dependent, are considered here constant 
for simplicity, i.e. u w u{k — kp) and w w w(k = kp), 

being the Fermi momentum. 
The electrons coming from the left lead acquire the time- 
dependent spin sensitive phase = n(^AB + o'^R{t)) 
(a = ±1), where &ab is the Aharonov-Bohm phase, 



while is the Aharonov-Casher phase produced by the 
modulation of the Rashba spin-orbit interaction on the 
ring = $° B , + sin(wt) 9]. The Hamiltonian given 

in Eq[T] can be rewritten by means of a plane wave ex- 
pansion in the following form: 

H{t) = H + e w sin(wi + ip) ^ d\d a + 

+ ^4 u [cos($°)A(t) -<TBin($°)B(t)]4cfcri +H.(g) 

ka 

where the static part Hq of the Hamiltonian is the same 
as in Eq.([T]) with: e{t) — > e , u — > uJ (tt^), §n{t) — > 
<I>5j, while the functions A(t) and B(t) are given by (see 
APPENDIX © : 



A(t) 



5(f) 



J2n{^ b Rj cos(2nwf) 



(3) 



J 2 „_ 1 (^)sin((2n-1)^), (4) 



where J n (x) are the Bessel functions of first kind. In 
absence of an external voltage, the quantum transport 
of particles through the structure is only due to absorp- 
tion/emission of energy quanta tku associated to the ex- 
ternal time-dependent fields. 

To calculate the pumped current we employ the 
Keldysh Green's functions technique as formulated in 
Ref . [17| . The current in the left lead is given by Ii(t) = 

- e Tt( N i) = r H H >*l]) ( N l = Ek^l^i). It can 
be rewritten in terms of the retarded and lesser Green's 
functions GrJ^ (t, t{) and of the advanced/retarded self- 
energies S"^, (ii, t) of the quantum dot according to the 



following expression: 



+ G< CT ,(Mi)£^v(^)}- 



(5) 



(6) 



For a time-dependent problem the retarded, advanced 
end lesser (r,a,<) Green's functions and respective self- 
energies depend explicitly on two time variables instead 
of one. Thus, employing the two time Fourier transform 
(see APPENDIX [Bj the current {5]) can be rewritten as: 



)( 



— TR P {[ dE l dE 2 dE 3 A{E 3 -E 1 
h^ Kei J (27T)3 

x[G r atr ,(E 1 ,E a )E< tr , a (E 2) E a ) + 
+G< a ,(E 1 ,E 2 )^l a , a (E 2 ,E 3 )}}, (7) 

where the lesser Green's function G < (E%, E 2 ) of the dot 
is given by the Keldysh equation: 



G^ {E\ , E 2 



(2tt)2 



G' r ( J B 1 ,ei)S < (a,6)G Q (6,S 2 ), 

(8) 
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and the following relation between the retarded and 
advanced quantities can be used: l E a (E 1 ,E 2 ) — 
[E r (E 2 , -Ei)] , where S = G or S. In order to compute 
the current the knowledge of G r , £ r , is required. Be- 
low we are going to calculate them for the non-interacting 
dot and the weakly interacting one. In both cases, the 
wide band limit (WBL) will be employed for simplicity. 

A. Non-interacting case (U = 0) 

In the U — case the expression of the self-energies 
can be obtained exactly By calling w — 7 it (7 6 R), 
the retarded and lesser self-energies can be expressed 
in terms of the corresponding Green's functions of the 
leads, namely [g r ka {t-,t')] sp = -i5 sp 9{t - t l )er t£ k^ t ~ t ) and 

[SfcjMOW = « v /(eg)e- te2(t -*' ) > by the following re- 
lations: 

Kf(tut2)= E 7 2 l«l a ^ < (*i,*a)W + 

+ J2 4|M| 2 cos($ s (t 1 ))cos($ p (t 2 ))[.gfei < (ti,i 2 )] sp (9) 

where s,p are spin indices (|, J,), /(efe) is the Fermi func- 
tion, while Eu = s r k = Ek in absence of voltage bias. 
In WBL limit and making the substitution J2k a ~ * 
J dep a {e), we get: 

KpihM) = -i6 sp 5(t 1 -t 2 )T { 1 2 /2 + 2cos 2 (<f> s (t 1 ))} 

E<(ti,ia) « i(W(*i -i2)r [ 7 2 +4cos 2 ($ s (ti))],(10) 

where we introduced the quantities r° = 2tt p\u\ 2 and 
/(fi-ia) = / ^f(e)e- t£{tl - t2 K By defining the left and 
right transition rates I^(ti)/2 = 2r° cos 2 ($ s (ti)) and 
T r s (ti)/2 = r°7 2 /2, the retarded self-energy can be writ- 
ten as ^ r sp (h,t 2 ) = -i8 sp 5{h - t 2 ){T\{t x )/2 + rj(ti)/2} 
and thus the retarded Green's function of the quantum 
dot takes the following formi^: 

G r sp (t,t') = gl p (t,t>) x 

x ex p{-^/ ; d*i[ri(*o+r:(*i)]}. 

sIp^O = -i6 sp e(t-t')exp{-i J ^ tftie(ti)}.(ll) 

It can be computed exactly after having performed 
the following plane-wave expansion of the retarded self- 
energy: 

K P (h,h) = -iS sp S(t 1 -t 2 )T°[i+^ 2 /2 + 

+ J (27r*S)cos(2$2) + 

+ cos{2<^° s )M(t 1 ) - ssin(2$°)A/'(i 1 )](12) 

where the following auxiliary functions have been intro- 
duced: 

00 

M(t) = 2^2 J 2n {2Tr<^f l ) cos(2nujt) (13) 

n=l 



Af(t) = 2E- 7 2n-i(27r$%)sm((2rc-l)urf), (14) 

n=l 

and whose Fourier transform is given by: 
V r sp {E x ,E 2 ) = -7r^{i^-fi)e; + 

00 

+iQ s 2 J2 J2n(2ir^)S(E 1 - E 2 +2nu;r])- Q| X 

77— ±l,n— 1 
00 

x E »7Ja»-i(2?r$£)$(£i - £ 2 + (2 ? i - l)ur])}. 

7)=±l,ra=l 

(15) 

An analogous plane- wave expansion can be performed for 
the lesser self-energy T,f p (ti,t 2 ), whose Fourier transform 
is the following: 

E< (E X ,E 2 ) = 2m6 ap {f(E 1 )6(E 1 - E 2 )Q\ + 

OO 

+ E Q^2„(27r^) x 

x/(£'i + 2nujr])5(E 1 - E 2 + 2nuj-q) + 

OO 

+ E iQ3VJ2n-i(2n<t>%) X 

77— ±l,n— 1 

x/(£i + (2n - l)^)^ - £ 2 + (2ti - 1)^)116) 
where are given by: 

Qi = 2r°{ 7 2 /2 + 1 + J (2^) cos(2$°)} 
Q 2 = 2T° cos(2$°) 

C| = 2r°ssin(2$°). (17) 

The substitution of Eq.® and of G r {E u E 2 ) in Eq.© 
permits to determine the G < (Ei, E 2 ). 

The knowledge of the retarded and lesser Green's func- 
tion enables us to calculate the current generated by the 
pumping procedure in the form of a trigonometric series, 
i.e. If(t) = Ig + E,T = iK cos(nwi) + < sm(nwt)], allow- 
ing to recognize the d.c. component of the current. 

B. Weakly interacting case (U ~ 0) 

The weakly interacting limit (U w 0) can be studied 
by means of a self-consistent Hartree-Fock theory which 
is known to give suitable results when the Coulomb inter- 
action U is smal l 20 i 21 (i.e., U -C r°). In this framework, 
the energy of the electrons on the dot is modified by a 
spin dependent term related to the occupation number of 
the electron of opposite spin and thus the spin dependent 
energy becomes e a (t) = e{t) + U (n s (t)) (a = —a). The 
occupation number (n a (t)) is calculated self-consistently 
by means of the relation i(n a (t)) = G< a (t,t). The re- 
tarded Green's function of the dot is modified by the 
interaction according to the expression: 

G r sp (t,t') = G r sp (t,t')e^{-iU J dtiingih))}, (18) 
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where G r sp (t,t') represents the retarded Green's function 
derived in the non-interacting case. In order to determine 
the interacting Green's function, we can write (n s (ti)) 
as a trigonometric series of sin(nwt), cos(nujt) with un- 
known coefficients calculated in a self-consistent way, as 
explained below. 



III. THE SINGLE PHOTON APPROXIMATION 

Hereon we focus on the weak pumping limit and con- 
sider only single photon processes, i.e. involving emis- 
sion or absorption of a single energy quantum hu>. The 
weak pumping limit (i.e. the case in which a pure sm((p) 
behavior of the current is expected) is very important 
in experiments where the higher harmonics contribution 
seem to be negligible even though others anomalies oc- 
cur. Such anomalies in the current-phase relation will be 
discussed here later. 



A. Non-interacting case (U = 0) 

Within the single photon approximation the self ener- 
gies (TT5l) - (fli)]) can be approximated as: 



7r6 sp {iS(E 1 -E 2 )Ql 



£<(£i,£ 2 ) « 2TrS sp {if(E 1 )Ql - g^Ji(27r$£) x 

x J2 + wMEi -E 2 + ^){19) 

v =±i 

The approximation is valid for small values of 27r$^ so 
that the inequality Ji(27r$£,) > J„(27r$^), n > 1, is 
verified. Using the above approximated form of (fT!)]) the 
retarded Green's function of the quantum dot can be 
rewritten as follows: 

GI«(M') = -i6 sp 6(t-t')x 

x exv{-i[e -iQi/2)}(t-t')Xx 

x exp|— i[Al dti sin(w£i) + 

+ A 2 d*icos(o;ii)]} ) (20) 

where we introduced the coefficients: 

Al = ff a; cos(^)+iJi(27r^)Q| 

A 2 = £ w sin(^), (21) 

and the upper index (1) stands for the single-photon ap- 
proximation. Making a further expansion of the retarded 
Green's function for small e^/cu leads to the result: 

G r s f\t,t') - -i5 sp 0(t-t') x 



x exp{-i[e -iQ^/2)](t-i')}x 

x [A s +AlC(t,t') + A s 2 S(t,t')], (22) 

where C(t,t') = cos(cji') — cos(ujt), S(t,t') = sin(wt') — 
s'm(ujt), while the coefficients A| have been defined as 
follows: 

Ag = Jo( f^M ) 2 J ( g^M ) 2 x 



x I ( 



A" 



(23) 



where I n (x) (n = 0,1) represents the modified Bessel 
function of first kind and order n. The above result 
can be conveniently rewritten in terms of the two-time 
Fourier transform and thus we have: 

G s y(E u E 2 ) - 2vrd sp j v^EA 

where we defined K* = (Af - ir]A a 2 )/2 and £> s (£:i) = 
-Bi — Eq + iQf/2. The knowledge of the retarded Green's 
function allows us to write the lesser Green's function by 
means of the Keldysh equation in this way: 

G< 1 HE 1 ,E 2 ) = 2iriQ s 1 ?° s (E 1 ,E 2 ) + 

+ 2^Q^J 1 (2^) vK n (Ei,E 2 ), 



T]=±l 

where we defined the following integral function: 



(25) 



^(E U E 2 



E 



(2tt) 2 as 



(26) 



The above function, disregarding terms quadratic in 
J\(x) and I\{x) and additional terms describing higher 
order processes (roughly cubic in ~ [V s (Ei)]^ 1 ), can be 
written in the simple form: 

~, (w r) _ f(E 1 +r ] u;)[Atf6(E 1 -E 2 +T 1 u) 



5 



Thus, the d.c. current generated by the time-varying 
parameters has the following final expression: 



(If (1) (*)> = 



eeiejAg[l-AS] f f(E)dE 
\V S (E)\ 2 
2eQ s 3 J 1 (2^%) 



LI •> 



J) = ±l ' 



un s 11 f(E)dE 



}.(28) 



Here is a coefficient obtained setting 7 = in Qf (for 
the left lead). The current (|28|) contains terms propor- 



where£>o(i?) = E — eo + iQl/2. By comparing the above 

expression with the G^jp (t, i) obtained by using (|30| one 
gets the following set of self-consistency equations: 



Q£ [ HE)dE_ 
2n J \Vg(E)\ 2 



G? 



,(2) - 



0. 



(33) 



tional to e 2 and ($'r) that can be interpreted as rec- induced shift- 



Once the above equations are solved, the d.c. current can 
be written as in Eq. (j2"8)) with the following interaction- 



tification terms, and terms proportional to £u^f> which 
contain information on the non-adiabatic pumping pro- 
cess, as will be clear below. 



e + Ua { °\ Af 



Af + Ua 



(i) 



C. The spin and charge currents 



B. Weakly interacting case (U ~ 0) 

To perform the analysis in the weakly interacting case, 
we consider as negligible the terms proportional to Uw. 
In addition we also consider U/lo as a small quantity, be- 
ing U of the same order of £ w , <&^. Within the Hartree- 
Fock theory we need to determine the energy of the quan- 
tum dot e a (t) = e + Un s (t) with n 9 {t) = {n s {t)). By 
using the single photon approximation, we write the oc- 
cupation number as a trigonometric series: 

n a {t) = 4°) + a^> sin(wi) + a {2) cos(wf), (29) 

(i) 

where the unknown coefficients a a have to be deter- 
mined self-consistently. In the interacting case the re- 
tarded Green's function takes the following form: 

G r £Xt,t') = -iS sp 9(t-t') x 

x exp{^[£ + CA4 0) -iQt/2)](t-t')}x 

x exp|— ^1 J dti sin(ujti) + 

+ Ai J dtxcos(ujti)]j, (30) 
where the coefficients Af have been defined as follows: 



\{ = £ w cos(v3)+zJ 1 (27r^)Q^ + [/4 1 
X'2 = e^, sin(<^) + Ua 



(2) 



(31) 



Note that since the coefficients ai appear as a factor 
of the interaction U , we have to calculate them only up 
to the zero order approximation in U and w. From the 
lesser Green's function obtained by (f30|) . we can write 
the occupation number in the following form: 



r(t) 



(A? 



<r\2 



[Qi - 2QlJ l {2-K^ R )sm{ujt)] 



2tt 
f(E)dE 

mm 2 



(32) 



To obtain an analytical expression of the d.c. cur- 
rent pumped in the left lead in the presence of a weak 
interaction and zero temperature, we need an approxi- 
mated expression for H* in (|28| in the limit of <C 1, 
and 271-$^ <C 1. In this limit one can use J\(x) ~ x/2, 
h(x) ~ x/2, J (x) ~ 1, I (x) ~ 1 (see APPENDIX [CJ 
and thus the coefficients take the following simplified 
form: 

V s „ kQrQI , suite) _ (e u cosjcp) + Ua { V) 

2lo + 1 2lo 2lo A ' 

while the quantity Ag[l — (Ag)] can be written as follows: 



Ag[l 



1 re 2 



a^/Je^ cos(<p) 



+0(1/lu 4 ). (35) 



Plugging these expressions in (|28p and performing the 
integral over the frequency, after an expansion up to the 
second order in u>, we can write the d.c. current in the 
single photon approximation (in units of 2T e/Ti) as fol- 
lows: 



2uj 2 



(2tt$£<£) 2 + 2ai ) Ue UJ cos(^) 



/i — e + t/(a, 



(o) „(oh 



°(u)\ 2 \\V° 



WW 
oj sin(</j) 



+ 



2|X>(//)| 



(m)I : 



-[(/*■ 



t/a. 



(0)n2 



2<£ cos((p)}, 



(?n 2 i 

(36) 



where the energies are measured in units of T°, while 
we defined q? = Qf /(2T°). The non-dimensional charge 
and spin currents, i.e. I c and I s , can be defined as 
I c = and I s — ^2 a <yia- The main feature of 

the expression for the d.c. current is the presence of 
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a non-sinusoidal current-phase relation already in weak- 
pumping. Indeed, contrarily to the adiabatic case char- 
acterized by a current-phase relation with definite odd 
parity (i.e. I c (—ip) = —I c (ip) ) in the time-dependent 
case any parity with respect to the sign reversal of ip is 
expected in the pumped current. This behavior is mainly 
related to finite frequency effects as well as to interaction 
effects. Eg. (1551) represents the main result of this work. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In order to make a comparison with the available ex- 
perimental data, we set T° ~ 10/ueV[22|. Such quantity 
is related to the dwell time Td by the following relation 
E T = h/Td ~ 2r°. Such quantity is relevant to define 
the various transport regimes at varying frequency uj. 
Indeed, for value of uiTd -C 1 one deals with the adiabatic 
regime, while in the opposite limit, i.e. ujTd 1, 
the non-adiabatic regime is approached. For typical 
experimental frequencies ranging from 10 MHz up to 20 
GHz, u!Td varies from ~ 10~ 2 up to order 10 and thus 
the MHz range of frequency can be safely considered as 
adiabatic. The adimensional frequency u which appears 
in Eq. (|36p is defined as ui = huj/T° — ujTd/ir. In this 
way a frequency of 25 MHz corresponds to uj — 0.01, 100 
MHz to uj = 0.04, 1 GHz to u = 0.4. In the following 
we study the behavior of charge and spin currents in 
the range of frequency uj £ [0.1, 0.5], thus our analysis is 
valid from adiabatic up to the moderate non-adiabatic 
limit. We also set the chemical potential fi as the zero 
of energy. From the analysis of the current i a , we notice 
the presence of two classes of terms contributing to 
the currents: 1) Terms proportional to e 2 or ($#) ; 
2) terms proportional to $^£,j. The first type of 
terms are non-adiabatic in nature. The second class of 
terms contains a term proportional to to, which can be 
recognized as the quantum pumping contribution, and 
a frequency independent term proportional to cos(ip) 
which can be interpreted as a rectification contribution. 
Such term is responsible for the non-sinusoidal behavior 
that leads to an anomalous current-phase relation as 
observed in Ref.[22j (page 3, first column, line 2). Very 
interestingly, the interaction effects also lead to a cosine 
term which is proportional to Ue^ R cos{p) /uj 2 . Such a 
term produces a deviation from the sinusoidal behavior 
also for small values of the energy U . Finally, the current 
%/j vanishes when the amplitude of the modulation e w 
and $^ go simultaneously to zero. 

In Figj2]the charge (dashed-dotted line) and spin (full 
line) currents, namely I c and I s , as a function of the 
phase difference ip between the time-varying parame- 
ters are reported for the following choice of parameters: 
7 = 0.05, $ AB = 0.49, <f>° R = 0.02, = 0.01, e = 0, 
e u = 0.025, u> = 0.1 and U — 0. A sinusoidal- like be- 
havior is observed even though the charge pumped for 
p = is different from zero and of the order 10~ 3 . This 



is a fingerprint of the anomalous current-phase relation, 
as discussed above. To put in evidence the dependence 




FIG. 2: Currents of charge (dashed-dotted line) and spin (full 
line) as a function of ip obtained for the following choice of 
parameters: 7 = 0.05, *ab = 0.49, $° R = 0.02, = 0.01, 
£0 = 0, e w = 0.025, ui = 0.1 and U = 0. 

on the interaction U, we present in Fig[3]the charge cur- 
rent computed at p — (dashed line) and p = 7r/2 (full 
line) as a function of U taking the remaining parame- 
ters as in Figf21 Smaller values of the interaction favours 
deviation from the sinusoidal behavior. Below we con- 
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FIG. 3: Currents of charge computed at ip — 0, I c (<p = 0) 
(dashed line), and tp — it/2, I c (<p ~ tt/2) (full line), as a 
function of U obtained for the following choice of parameters: 
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0.05, $ab = 0.49, 
0.025, u = 0.1. 



0.02, $ 



R 



0.01, £ = 0, 



centrate on the role of spin-orbit interaction and choose 
the Aharonov-Bohm flux close to half integer values in 
unit of the flux quantum $0 = h/e where the charge 
current is activated by photon-assisted tunneling (PAT) . 
Away from the above values of the Aharonov-Bohm flux 
the currents present an oscillating behavior as a function 
of the applied magnetic flux $ab similar to the one al- 
ready discussed in a previous work^ 3 -. 
In Figf4]we plot I c (dashed-dotted line) and I s (full line) 
as a function of the static Aharonov-Casher phase Q a R for 
pumping frequency uj — 0.2, ui — 0.3, uj = 0.4 (from top 
to bottom) and by setting the remaining parameters as 
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FIG. 4: Charge (dashed-dotted line) and spin (full line) cur- 
rents as a function of <E>5j obtained for the following choice of 
parameters: 7 = 0.05, $ab = 0.49, = 0.01, e = -0.025, 
Euj = 0.05, tp — 57r/4, U — 0. The upper panel is obtained for 
u> = 0.2, the middle panel for uj — 0.3 and the lower panel for 
uj = 0.4. 

frequency from uj = 0.2 (500 MHz) up to 0.4 (1 GHz) 
zeros of the charge currents start to appear and thus it is 
possible to obtain pure spin currents in the non-adiabatic 
regime. It is worth mentioning that currents of ampli- 
tude 10 -2 in dimensionless units correspond to ~ 50pA 
in dimensional unit, thus the pure spin current we find is 
sizable. 



To analyze the role of a weak Coulomb interaction, we 
plot in the upper panel of Fig [5] the charge and spin cur- 
rents for the same parameters as in FigH] (to = 0.3) and 
by setting U = 0.02. A qualitatively different behavior of 
the charge current as a function of the Aharonov-Casher 
flux is observed. In particular, when the interaction en- 




0.04 0.06 




FIG. 5: Charge (dashed-dotted line) and spin (full line) cur- 
rents as a function of obtained for the following choice of 
parameters: 7 = 0.05, = 0.01, eo = -0.025, = 0.05, 
uj = 5tt/4, uj = 0.3, U = 0.02 and $as = 0.49 in the upper 
panel, $ab = 0.52 in the lower panel. 

ergy U is of the same order of magnitude of the pumping 
frequency uj, additional zeros of the charge current ap- 
pear and this is a very appealing situation for spintron- 
ics devices. For instance, looking at Figf5] (upper panel), 
one observes a pure spin current close to = 0.015 and 
0.03. 

In the lower panel of Figf5]we plot charge and spin cur- 
rents as done in the upper panel and by setting the 
Aharonov-Bohm flux to §ab — 0.52. In this char- 
acteristic oscillating behavior of the currents controlled 
by using a magnetic flux is visible. 

Another interesting phenomenon is the asymmetric 
contribution to the current of the photon absorption and 
emission as a function of the dot level eo 1 as also reported 
in Ref . [24] . When the dot level lies above the Fermi level 
(eq > 0), an electron on the dot can jump in the left 
lead by emitting a photon. For £0 < 0, an electron on 
the dot can reach the left lead only by means of the ab- 
sorption of a photon since no voltage bias or temperature 
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gradient are present. Because of the interference between 
these two-photon sources, boson-assisted tunneling onto 
the dot gets suppressed while tunneling out of the quan- 
tum dot is enhanced. The asymmetric behavior of the 
d.c. current as a function of the dot level is shown in 
FigEl In the upper panel, we plot the charge current I c 



the upper panel of Fig(7] we plot the current obtained 
for e u — 0.05, while in the lower panel this parameter 
is set to zero (single-parameter pump). By comparing 
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FIG. 6: Charge current I c as a function of £o obtained for 
the following choice of parameters: 7 = 0.05, <&ab = 0.49, 
<f>° R = 0.05, = 0.01, if = tt/2, U = and e w = 0.05 for 
upper panel, g w = for lower panel. Each panel contains 
curves obtained for u) = 0.1 (dashed line), ui = 0.25 (dashed- 
dotted line), u) = 0.5 (full line). 

as a function of the dot level £0 and by setting the remain- 
ing parameters as: 7 = 0.05, <&ab = 0.49, <S>° R = 0.05, 
= 0.01, <p = n/2, U = and e u = 0.05. As can be 
seen, when the frequency lo is increased from 0.1 (dashed 
line) up to 0.5 (full line) a strong peak is formed at Fermi 
energy and the asymmetry of the current with respect to 
the £o = becomes more evident. It is worth to mention 
that, since the relative phase cp is 7r/2, all the terms in 
the current proportional to cos(<p) are suppressed, while 
the pumping term takes its maximum value. In the lower 
panel we set e u = 0, while the remaining parameters are 
fixed as in the upper panel. In this case the device works 
as a single-parameter pump associated to the Aharonov- 
Casher flux and the current is proportional to (Q R ) 2 . A 
comparison between the upper and the lower panel shows 
that the pumping mechanism is the dominant one at the 
Fermi energy. Furthermore, we have verified that a small 
interaction does not alter too much the picture given so 
far. 

The same analysis performed in Figs[5] can be repeated 
by setting ip = to include the cos(</?) contribution. In 
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FIG. 7: Charge current I c as a function of £0 obtained for 
the following choice of parameters: 7 = 0.05, §ab = 0.49, 
$° R = 0.05, = 0.01, Lp = 0, U = and e w = 0.05 for upper 
panel, — for lower panel. Each panel contains curves 
obtained for lo = 0.1 (dashed line), lu — 0.25 (dashed-dotted 
line), lo = 0.5 (full line). 



the results, one observes an enhancement of the abso- 
lute value of the high frequency currents in the case of 
double-parameter modulation (upper panel) and close to 
the Fermi energy. 

From the analysis above one observes that, within the 
considered parameters region, the dominant mechanism 
for the generation of the d.c. current is the finite fre- 
quency quantum pumping. Indeed, close to the Fermi en- 
ergy such currents take values which range from ~ 70pA 
up to ~ 190pA (see the upper panel of Figs[6]), while 
in the other cases the generated currents present values 
of about 10% of those induced by the pumping process. 
Thus, for Q ab close to half-integer values the quantum 
pumping induces the main contribution to the current, 
while away from this flux region the rectification currents 
are dominant. 



V. CONCLUSIONS 

We studied the time-dependent charge and spin trans- 
port (pumping) in a Aharonov-Bohm-Casher ring se- 
quentially coupled to a weakly interacting quantum dot 
by using a non-equilibrium Green's function approach. 
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By varying a considerable number of parameters, we 
showed that the proposed device can work as a spin cur- 
rent generator and analyzed all its characteristics, in- 
cluding rectification effects. When the energy level e(i) 
on the dot and the Aharonov-Casher flux are periodically 
modulated in time with a frequency u>, a d.c. current is 
observed in the leads. Contrarily to the adiabatic case, 
the current-phase relation presents two additional cosine 
terms: The first one comes from the interaction on the 
dot, while the second can be interpreted as a rectifica- 
tion effect, as already noted in Ref.[22|]. We also showed 
that Coulomb interaction effects can enhance the rectifi- 
cation effects. As a function of the spin-orbit interaction 
and close to the non-adiabatic regime, the results of the 
charge current show the appearance of additional zeros 
at varying the Aharonov-Casher flux. Thus, the finite 
frequency regime close to 750 MHz (ui = 0.3) is suitable 
to obtain pure spin currents useful in spintronics. Such 
currents are of the order of magnitude of ~ lOOpA as 
detected in the experiments on quantum dots^ 5 .. Finally, 
the analysis as a function of the dot level showed a char- 
acteristic asymmetric behavior and the comparison be- 
tween the single parameter pump and double-parameters 
one showed a considerable increase of the d.c. current 
in the second case. The proposed device can be eas- 
ily fabricated on a GaAs/AlGaAs two-dimensional elec- 
tron gas using e-beam lithography to define the ring and 
dot region modifying, for instance, the system studied in 
Ref.fH. 



cos(a;sin(0)) = Jp(x) + 2 J 2n (x) cos(2?- 



n=l 



exp{Acos(<9)} = 4(A)exp(m0). (Al) 
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APPENDIX B: TWO-TIME FOURIER 
TRANSFORM 



The two-time Fourier transform has been defined ac- 
cording to the following definitions: 



dEi dE 2 

x 

2tt 2tt 

x g(E 1 ,E 2 )exp{-iE 1 t 1 + iE 2 t 2 } 



g(E 1 ,E 2 ) = f dhdhgfahjexpiiEttt -iE 2 t 2 }Bl) 



APPENDIX C: APPROXIMATION OF J {x), J^x), 
Io(x), Ii(x) FOR x w 0. 
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APPENDIX A: BESSEL EXPANSION 

Throughout the paper the following expansions have 
been exploited: 



Mx) =< l- — + 0(x 3 ) 

j l{x) ~ ?L + 0(x 3 ) 

Io(x) ~ l + ^ + 0(x 3 ) 

I x {x) ~ | + 0(x 3 ). 



(CI) 



sin(:rsm(0)) = 2^ J 2n _i(a;) sin((2n - 1)0) 
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